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Abstract 


A perturbative QCD calculation of heavy flavor quark fragmentation into 
heavy flavor baryons is developed along the lines of corresponding heavy me¬ 
son models. The non-perturbative formation of the baryon is accomplished 
by implementing the quark-diquark model of the baryons. Diquark color form 
factors are used to enable the integration over the virtual heavy quark mo¬ 
mentum. The resulting spin independent functions for charmed and bottom 
quarks to fragment into charmed and bottom baryons with spin 1/2 and 3/2 
are compared with recent data. Predictions are made for the spin dependent 
fragmentation functions as well, particularly for the functions gi and hi in 
the case of spin 1/2 baryons. 
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I. INTRODUCTION 


Quarks produced in high energy processes materialize by evolving into jets of hadrons. 
The particular hadronic fragments and their kinematic dependences are of considerable in¬ 
terest. This fragmentation process reveals the features of the non-perturbative regime of 
QCD. The inclusive process of quark fragmentation into a single observed hadron along 
with any number of unobserved accompanying particles is described by a set of fragmen¬ 
tation functions. These fragmentation functions, defined in terms of appropriate light cone 
variables, kinematic variables and invariants, are probability distributions. These functions 
have received considerable attention in recent years. While experimental information be¬ 
yond the pion distribution (presumably from light quarks) has been slow in accumulating, 
theoretical interest has been growing. The particular functional form for heavy flavored 
quarks to fragment into heavy flavored hadrons is of special interest. 

Experimentally it is possible to identify heavy flavored jets by the production and char¬ 
acteristic flavor changing weak decays of the hadrons. Theoretically this situation has been 
studied using Operator Product Expansion techniques, light cone quantization, QCD per¬ 
turbation theory, and Heavy Quark Effective Theory, among other methods. Some of these 
methods yield general properties that reflect the overall structure of QCD, as it is currently 
understood. Other approaches take particular models of the low energy behavior expected 
from QCD, but in regions that are not perturbatively calculable. The particulars of the var¬ 
ious approaches are near the point of being tested against experiment. One general feature 
is known - the peak of the hadron distribution moves toward higher momenta as the quark 
mass increases. This feature is a result of the kinematics implicit in most models of the 
non-perturbative process, and is incorporated in the phenomenological Peterson function 
that is used by experimenters to £t the sparse data on heavy quark fragmentation |^. 

It is more difficult to test the spin dependences of the fragmentation processes experi¬ 
mentally. Yet these dependences are very important to know. They reflect the details of 
the primarily non-perturbative mechanism by which parton polarization is passed on to the 
hadrons. The spin-dependent fragmentation involves the reverse of the process by which the 
nucleon spin is shared by its partons (the“spin crisis”), and may reveal a similarly mysterious 
decoupling of valence quark spin and hadron spin for some regions of kinematics. 

For the fragmentation of a heavy flavor quark into “doubly heavy” mesons, perturbative 
QCD may provide a starting point to a theoretical determination of the fragmentation 
functions. As examples, the fragmentation of a c-quark into the J/T or the b-quark into 
the Be, were calculated several years ago The reasoning follows the observation that 

a heavy quark or diquark pair must be produced to form the hnal hadron. At least one 
member of the pair will be nearly on-shell because the heavy flavor hadrons have small 
binding energies relative to their masses. So, the gluon producing the heavy pair must 
carry large squared time-like 4-momentum, This gluon will be shaken off by the virtual 
fragmenting quark and the relevent coupling will be of order as{k‘^), which will be small. 
Hence perturbative QCD will be applicable. Non-perturbative effects will be incorporated 
into the binding of the initial quark with the pair produced heavy quark or diquark. If 
this is the case, the fragmentation functions are calculable, at an appropriate scale. The 
parton shower that accompanies the jet is a result of QCD radiative corrections, which can 
be obtained from the Renormalization Group or the Altarelli-Parisi equations. 
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Such calculations have been performed and scrutinized. It has been shown that in the 
heavy quark limit (i.e. the mass goes to inhnity) the functions have the form expected 
from more general considerations |^. This corresponds to the heavy meson taking all of the 
heavy quark’s momentum; the distribution becomes a delta function at z = 1. The l/mg 
corrections are calculated also. In any case, this approach can predict the spin-dependent 
fragmentation functions along with their momentum and mass dependences. In the heavy 
mass limit, of course, the spin of the heavy quark is conserved, so the spin dependence is 
simple. What is of phenomenological interest is the next order correction, at least, since 
that has non-trivial spin dependence. 

The spin dependences of fragmentation are most readily studied experimentally by ob¬ 
serving baryons rather than mesons. This is true for the production of hyperons or heavy 
hyperons (Ac, A;,, etc.), wherein the weak, parity violating decays provide polarization analy¬ 
ses P]. To consider fragmentation into baryons in this perturbative scheme, the three quark 
system has to be confronted. A simple alternative is to consider the baryons as quark- 
diquark bound states 0, and to use the same perturbative method as for the mesons. In 
order for the perturbative calculation to be useful the creation of a heavy pair of quarks or 
diquarks must be an intermediate step. Ideally then, doubly heavy baryon fragmentation 
would be an appropriate testing ground for this scheme. Such data is sparse, however. 

To begin to see the structure it will be worthwhile to stretch the region of applicability to 
the “singly” heavy baryons. We have been carrying out this program to see the expected spin 
and kinematic dependences, with the hope of providing an experimentally testable model [§ . 
A similar approach has been developed independently by Martynenko and Saleev [^, but 
with signihcantly different assumptions about how to represent the diquark structure. Of 
immediate interest is the question of whether the baryon fragmentation functions have the 
same kinematic dependence as the meson case. In general the answer is no in this model. 
Secondly, does this approach give the right magnitude for fragmentation into heavy flavor 
baryons? With our careful specihcation of the diquark chromodynamic form factors (unlike 
Martynenko and Saleev) the answer is yes. The spin dependent fragmentation is interestingly 
distinct from the naive heavy quark limit in detail. The calculations and results will be 
presented below, along with a comparison with some recent data. 


II. PERTURBATIVE CALCULATION 

The hrst calculations of the fragmentation functions in the perturbative scheme were ap¬ 
plied to some of the inclusive heavy flavor meson decays of the as produced at LEP P,||. 
The partial width for the inclusive decay process ^ H + X can be written in general for 
any hadron H as 

dr(Z° ^ H{E) + AT) = V I'dz dt{Z^ ^ t{E/z) + X, p) fJ.), (1) 

where H is the hadron of energy E and longitudinal momentum fraction z relative to the 
parton i, while /i is the arbitrary scale whose value will be chosen to avoid large logarithms. 
The quark and the hadron can carry spin labels as well, and appropriate spin-dependent 
fragmentation functions will be included, as we will show later. The fragmentation function 
fi) enters here in a factorized form (that can be maintained through the evolution 
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equations). Upon obtaining the in the model to be described, its evolution to 

observable scales is developed through summing leading logarithms via evolution equations. 

Now, consider the hnal state with one heavy flavor meson, say the Be for definiteness. 
We will soon replace this process by one involving a heavy flavor baryon. To leading order 
the Be meson arises from the production of a pair of b-quarks, in which one of the quarks 
fragments into the meson. As Fig. illustrates, with Q = b,Q' = c,Q' = c, the perturbative 
contribution involves the virtual b-quark radiating a hard gluon (We work in axial gauge, 
so that there is no contribution from the opposite quark). The hard gluon produces a heavy 
flavor pair of c-quarks. The gluon must have energy at least twice the charm mass, since 
the c and c are both on or near their mass shell. So the coupling as{k‘^) is small, justifying 
the perturbative approach. 

For nearly matching 4-velocities {v = p/m) the b and c form the Be meson bound 
state at roughly the same 4-velocity, with amplitude given by various projection operators 
multiplying the Bethe-Salpeter wave function x(p, qr) {qr is the p{b)—p{c) relative momentum 
while p is their sum). Since the doubly heavy mesons are weakly bound objects (the sum 
of constituents’ masses are near the bound state mass), the wave function is expected to 
dampen non-zero relative 3-momenta qr in the hadron rest frame [^, so that the b and 
c 4-velocities are fixed at p/M. Then the integration over the relative momentum of the 
two heavy quarks in the full decay probability can be replaced by the squared amplitude 
evaluated at equal 4-velocities (for the two constituents and the hadron). The remaining 
integration over qr applies to \x{p, <?r)P which yields the square of the wave function at the 
origin in the hadron rest frame. That wave function is known from non-relativistic quark 
models for the heavy-heavy meson system, or, more directly, from the meson-to-vacuum 
decay constant. 

The same procedure can be applied directly to the baryons, if the quark-diquark model 
of the baryons is used. The hard gluon in the process must produce a diquark-anti-diquark 
pair, D — D, and the diquark (color anti-triplet) combines with the heavy flavor quark Q to 
form the baryon Bq. The relevent wave function will be calculated from a fairly successful 
quark-diquark model 0. 

Note that an alternative scenario has been proposed by Falk, et ai, in which the heavy 
quark fragments into a heavy diquark first, and then the diquark dresses itself to form 
the baryon 0 with probability of one for the latter. This leads to very different results, 
as pointed out in Ref. 0. This scenario will not be used here, since the processes we 
are studying involve diquarks that do not necessarily carry the heavy flavor of the quark. 
The latter authors PJT^ have performed calculations that are similar in spirit to part of the 
procedure we follow below, although not emphasizing the spin dependent structure functions 
that we calculate below. 

The tree level amplitude for Fig. |l], Ai, can be evaluated explicitly from perturbation 
theory. The decay rate for unpolarized ^ Be + c + b ot Bq + Q + D, each an exclusive 
channel, can be written generically as 


Fi = 


1 


2Mp 


'] (27r) V(Z -q-p- p')-J2\Ai 


( 2 ) 


where q, p, and p' are the 4-momenta of the b, Be and c (or the Q, Bq and D), respectively, 
and I Alp is summed and averaged over unobserved spins and colors. We use the notation 
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[dp] = d^p/{IQ tt^po) for the invariant phase space element. In spin dependent fragmentation 
the sum will only cover the unobserved outgoing parton spin labels (c or D in this explicit 
case). To isolate the fragmentation function, the production of the fragmenting quark (dh 
of Eqn. ID) must be factored out. The hctitious decay width for the ^ h + h oi Q + Q, 
with the b or Q-quark on shell is 

To = ^^J[dq][dq]{27r)^6\Z (3) 

with q the b or Q-quark 4-momentum. 

To obtain the full inclusive width the unobserved quark degrees of freedom must be 
integrated over. By introducing the variables q, the off-shell quark’s 4-momentum, and 
s = q"^, the square of the virtual mass, the two body phase space for p and p' can be written 
as an integration over ^ = {po + Pl)/( qo + Ql) and s. Note that the transverse momentum of 
the hadron, px, and the unobserved quark, p'x = — Pt (relative to the fragmenting quark 
momentum), are hxed for each pair of z and s values via the relation 


s = q^ = {qo + q3){qo - qs) 


+ Pt^ , Pt^ 

i ^ 1-z 


( 4 ) 


where M and m' are the masses of the hadron and the unobserved quark or anti-diquark, 
respectively. Then the phase space integration in Eq. ^ can be written for the on-shell 
hadron and unobserved quark and anti-diquark production as follows: 


j[dq\[dp\[dp] {2 t:Y 5^{Z -q-p-p') 

7r)V(Z -q-q) j[dp][dp'] (27r)^(5^(g - p - p') 

= JdsJ[dq][dq]{27iY6^{Z - q - q) 

The variables po and ps have been replaced by s and z, and the integration over px has been 
performed via the delta function that requires 

Pt‘^ = z{l - z)[s - sth], ( 6 ) 

where s+/, = — -t- is the minimum value that s can assume. Note that the azimuthal 
integration in px has been performed assuming there is no such dependence in the amplitude. 
This will be true for spin averaged probabilities and for products of helicity amplitudes, but 
not for other orientations of quark or hadron spin. If spin projection operators are used in 
trace expressions, care must be taken to integrate out the azimuthal dependence Erst. The 
integrations over q and q will be common to the direct production of an on-shell quark in Eq 
and the off-shell quark that fragments in Fi. Providing the production dependence (l^oP) 
can be factored out of the full probability (jAip), this will allow the fragmentation process 
to be defined irrespective of the production mechanism, obviously an essential feature of any 
model. The factorization will be possible in the appropriate large momentum limit. 

To match the integrand to the fragmentation function of Eq |I| the integration will be 
performed over the variable s, keeping ^ fixed and letting Mz and q^ —>■ oo. The s integration 
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ranges from Sth to {Mz — rng)"^. Since the gluon propagator in the amplitude emphasizes 
low values of and the heavy quark propagator favors s not far from on-shell; the major 
contribution to the s integration appears at low s. Thence, the upper limit of the integration 
over s can be taken to cx) to facilitate the evaluation of the dehnite integral. In the large 
Mz or go oo approximation the transverse momentum of the hadron is small relative 
to po and ps, since Eq. ^ shows the transverse momentum is independent of go at hxed s 
and 2 ;. Thus it is sensible to ignore the transverse momentum in the relation p = zq (after 
carefully evaluating s dependent terms in the integrand). Once the square of the amplitude 
Ai is summed over spins and simplihed by dropping non-leading contributions, the width 
for > QQ can be factored out of the expression Eq. via 


Dq^h{z) 


-T lim 



1^1 


1^0 


( 7 ) 


leaving an integral over the fragmentation function, since the production probability for the 
relevent quark has been factored out. Then 


Dq^h{z) 


27 nmq 


dsF{z, s), 


'Sth 


( 8 ) 


where i?(0) is the Bethe-Salpeter wavefunction at the origin, and F{z,s) is the remaining 
integrand, which depends on s = g^, 2 ; and the quark masses, with the Qq —>■ 00 having 
been implemented. So the partial width for —y H + X is given by an integral over the 
virtuality of the heavy quark and the phase space of the unobserved degrees of freedom. 

We now proceed with the calculation of fragmentation functions for (singly) heavy flavor 
baryons. The basic covariant coupling of diquarks to gluons was written long ago |0]. There 
is one coupling constant for the scalar diquark color octet vector current coupling to the 
gluon held—a color charge strength, along with a possible form factor Fg. The momentum 
space color octet current (which couples to the gluon held vector) is 

= 9.F.(P)(p + pXS^iXi^S", (9) 


where p and p' are the scalar diquark 4-momenta and k = p' — p. For the vector di¬ 
quark there are three constants - color charge, anomalous chromomagnetic dipole moment 
K, and chromoelectric quadrupole moment A, along with the corresponding form factors, 
Fb, Fm, andFg. 

JX' = ■ f'’'(p')l(p + p% (10) 

+(1 + K)FM(«:^)|e“(p)p • ("Hp') + tl'(p')p' ■ e“(p)| (11) 

+AFQ[k%t'i(p)4'(p') + ivi“(P) ■ /*(p')]*:''V(p + p\ } , (12) 

where A is the color octet index, a, (3, ..., are color anti-triplet indices, the e’s are polarization 
4-vectors for the diquarks. 

In the perturbative diagrams involved here, the virtual heavy quark emits a time-like 
oh-shell gluon, that,in turn, produces a diquark-antidiquark pair while attaining nearly on- 
shell 4-momentum. The diquark combines with the heavy quark to form a heavy flavor 
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baryon, whose amplitude for formation is related to the Bethe-Salpeter wavefunction for 
the diquark-quark system. As in the meson production calculations, it is assumed that the 
constituents are heavy enough so that the binding is relatively weak, i.e. the quark and 
diquark are both on-shell and the binding energy is negligibly small. This is expected to 
be true for constituents with masses well above Aqcd, and even the light flavor diquarks 
almost satisfy this constraint. The basic perturbative amplitude is shown in Fig. |l| with the 
Q'-quark line replaced by an (anti-)diquark D line. 

It should be realized that the integration (over s, the square of the virtual heavy quark 
mass) involved in the calculation would diverge for point-like vector diquarks, since the gluon 
coupling to a pair, Eq. 13 carries momentum factors. The virtual mass in the integration, 
y/s, is passed on to the gluon and, subsequently, to the gluon-diquark vertex. Hence it 
is essential to regulate the integrand by some means. This is best accomplished via the 
chromoelectromagnetic form factors for the gluon coupling to the diquark. The form factor 
approach makes physical sense - it is a result of the compositeness of the diquarks. And for 
consistency, once the vector has form factors, the scalar diquark must have one also. 

There is no direct information about the chromoelectromagnetic form factors. We may 
expect that the ordinary electromagnetic form factors will have the same functional form as 
their QCD counterparts—the source of both sets of form factors is the matrix element of a 
conserved vector current operator. In the relevent case here, though, the vector operator is 
the gluon held — a color octet. Also, what is of concern here is the time-like region of the 
form factor. For diquarks, of course, there is not any direct empirical evidence about their 
electromagnetic form factors, but diquark-quark models of the nucleon have constrained the 
parameterization of the form factors. For one thing, the dimensional counting rules lead to 
l/|g|^ asymptotic behavior of nucleon form factors (at asymptotic momentum transfer the 
baryon is a three quark system). A quark-diquark nucleon must approach this asymptotic 
behavior also, for consisitency. For a point-like scalar diquark bound to a quark, the asymp¬ 
totic behavior will be l/|q'P from dimensional counting for the exclusive pair production. 
Hence the composite scalar diquark must have an effective coupling to the gluon, i.e. a form 
factor, that approaches asymptotia as l/Igp. Since the vector particle has a polarization 
4-vector associated with it, an extra power of momentum arises in the asymptotic amplitude. 
It becomes necessary for the charge and magnetic form factors to have 1/|5'|^ asymptotic 
dependence, and the quadrupole l/|g|® behavior. 

Using a quark-diquark model of the nucleon, and the well measured electromagnetic form 
factors of the nucleon, Kroll, et al. have obtained electromagnetic form factors for the 


diquarks. Two vector form factors are given 1/|5'|^ asymptotic behavior (the third is set to 
zero) and the scalar behaves as l/lijp. For the nucleon form factor study 0 as well as a 
recent study of higher twist contributions to the nucleon structure functions |]I^, the scalar 
diquark form factor and vector diquark form factor are assumed to have simple pole and 
dipole forms, respectively, with pole positions Ms and My above 1 GeV, 


2 \2 


Fs{k^) = 1/(1 - e/Ml), FE{k^) = 1/(1 - k^/M^) 

FM{k^) = il + n)FE{k^), FQ{k^) = 0. 


(13) 


The pole position values are somewhat higher than the dipole position for the overall nucleon 
form factors - near 800 MeV If we make the assumption that the color form factors have 
the same functional form as the electromagnetic form factors, we can proceed. However,the 
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region of most relevence for the fragmentation functions is time-like below the 4m^ 
threshold. In the s integration that will be performed here, the time-like region begins 
at ^m^Diquark value z = 1/(1 -|- mci/mB), and at higher values for other choices of 

2 ;. This implies that the integration region either overlaps or comes near to overlapping 
the pole positions. The pole singularities have to be tamed, and the hnal integration may 
be very dependent on the method used to moderate the singularities. Treating the poles 
as real resonance positions, including a small imaginary part, on the order of the nearby 
vector meson width, would be sensible physically. However, since the color octet form factor 
would be dominated by color octet vector mesons, and the latter are not expected to be 
strongly bound or narrow resonances, pole positions with large widths may be preferred. 
This would hide our ignorance and provide an interpolation between the space-like and 
time-like asymptotic regions. That is the ansatz we adopt. 

The neccessity for diquark form factors has an important consequence theoretically. In 
the light cone expansion, the baryon production via three quarks would contribute to the 
leading twist fragmentation functions. Dimensional counting requires the behavior to 

which we alluded above. But the diquarks depart from this behavior except at asymptotia. 
Hence the diquark form factors produce non-leading twist behavior for the fragmentation 
functions. Gluon contributions are buried in those form factors. In terms of the full set 
of such functions |^|, we have more than just /i, gi and hi. There are non-leading twist 
functions like ei and h 2 that will receive contributions at next-to-leading twist. When we 
determine what we call /i, gi and hi, we actually have some non-leading twist contributions 
that have not been disentangled. 

The amplitudes for the baryon production can now be calculated. The spin 1/2 ground 
state baryons are composed of a scalar diquark and a heavy quark in an s-state. There is 
only one coupling, and it involves the Fs- The amplitude is 


^Sl/2 — 


^/(O) 




Fs{,k‘^)UBgs[k\ - 2mdVx]P^ 


(14) 


where 


pA 




771^(1 -h V) -h k 


T. 


m 


Q) 


(15) 


For the vector diquark baryons, there are two form factors (we take the quadrupole to be 
zero - it falls as 1/|5'|® asymptotically). The chromomagnetic coupling involves a parameter 
K, the “anomalous chromomagnetic moment”. This is taken to be -1.10, as will be explained 
below in Section IV on comparing with data. The s-state baryons are spin 3/2 and 1/2, 
which we will refer to as 1/2'. The 1/2' lies between the 3/2 and the ground state 1/2 baryon. 
The amplitude for vector diquarks to be produced, along with the heavy quark, contributes 
to both 3/2 and 1/2' states. The amplitude is conveniently divided into a chromoelectric 
and chromomagnetic part, involving the two distinct form factors. The chromoelectric part 
contributing to the spin 1/2' baryon is 


Aei/2 = - [fcA - 2mdVx]P^. (16) 

The chromomagnetic contribution to the spin 1/2' baryon is 











( 17 ) 


Amiii = + k)(7b757'‘—- AK]P^- 

For the spin 3/2 baryon the corresponding amplitudes are 

Ae 3/2 = -^^^FE{k^)n9sel[kx - 2m,vx\P\ (18) 

and 


A 


M3/2 


^(0) 

x/2md 


FEik"^)^! + K)'^Bgs[gf,\{,e^v)md - e\k^j\P^ 


(19) 


Each amplitude should be multiplied by the color factor 4/3v^. In these amplitudes, ■^(0) 
is the Bethe-Salpeter wavefunction at the origin (for the s-state Q-diquark system), rud is 
the appropriate diquark mass, f/s is a spin 1/2 Dirac spinor for the baryon, e^{p') is the 
polarization 4-vector for the unobserved anti-diquark, is the Rarita-Schwinger spinor for 

the spin 3/2 baryon, v = p/M is the 4-velocity for the heavy baryon of mass M, F is the 
production vertex for the heavy quark-antiquark pair, k is the 4-momentum of the gluon 
and, with n = (1, 0, 0, —1), the corresponding propagator in axial gauge is 




1 / uu ^^k'^ + k^rP 

— [Q^ - 

{nk) 


( 20 ) 


Considerable simplihcation of these amplitudes follows. Recall that the formation of the 
baryon requires that the quark and diquark carry the same 4-velocity as the baryon, so that 
k = rp + p', with r = m£)/M and mq -|- rriE = M in the weak coupling approximation. This 
leads to many simplifying relations among the kinematic variables. Of particular importance 
is the relation = r(s—mg), which ties the gluon propagator to the heavy quark propagator 
as the s integration is performed. The resulting simplihed forms for each of the amplitudes 
become 
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^51/2 

AeI/2 




A-EZI2 


^M3/2 


V'(O) ^ ,, 2 ^^^ o 2 /2M2(1 - r) + M /; 


A/2m^ 

V’(O) 


Fs{k^)UB2gi 


[np) 


(s - mlY 


{nk){s - ml) 


r. 


( 21 ) 


V3mrf 


FE{k^)UB75 


2gl 


M{s - ml 


2 'i2 


I(£'p) + M ^'] 


(np) 


*[2M"(1 - r) - 2^{kp) + M /:]r, 
[nk) 


( 22 ) 


FE{t!‘){l + k)Ub1!, 


rM(s — ml 


2 'i2 


{—2(A;p)(pe'^)(l — r) 


V’(O) 

yj^rrid 

+2j^{kpf - 2r{pe^)p^{kp) + (3r - 2)M(pe^) /t + 2Mj^{kp) /: 

[TlfZ ] [Tin/ ) (TZfC) 

—2rM{pe^){kp) ^ + 2r{kp)M — (kp) /(jr, 

[TiKi ] 

^(0) c. J2M\1-t) + M 


(23) 


a/2^ 




V (s-"^q)^ 


(np) 


(nfc)(s-m^) 


(24) 


V’(o) 

V2md 


FEik^)ii + K)ngl 


r(s-ml)^ 


71 

{—2r{pe^){kp) ^ 


nk) 


-2(1 - r)k^(pi') + 2(kp)‘^^kE - / A}r. 


(25) 


There are 3 cases to consider for each heavy quark flavor—3 final state baryons. For the 
two states resulting from the vector diquark, the electric and magnetic amplitudes must be 
added together. Then for each baryon, the amplitude is squared and a trace is taken to sum 
over spins (including spin projection operators for the spin dependent cases). The analog of 
Eq. |] is obtained for each baryon. By carefully organizing the terms in the integrand, the 
width for the inclusive production of the virtual heavy quark can be divided out to yield 
the analog of Eq. for each baryon. Finally the integration over s = can be performed 
numerically—the form factors make it difficult to write an analytic expression for each case. 
The resulting z dependent fragmentation functions are the “boundary” functions, obtained 
at a scale p? at the threshold (2mTmg)^, with the strong coupling constant a{p = 2m]:,). 
To consider higher momentum scales, the Altarelli-Parisi evolution equations are used. This 
strategy essentially sums the leading log contribution of the parton shower that is generated 
by the off-shell heavy quark. 

We have taken some particular cases to illustrate the results. For the c(su) or c(sd) 
baryons, the Sc states, the diquark is given a mass of 0.9 GeV/c^ and the ratio of diquark 
to hadron mass is r = 0.33. Form factor parameterizations are discussed in Section 4. The 
resulting function, fi{z,Q‘^) is shown in Fig. || for the boundary value at p = 2mD + mg 
and for Qo = 5.5 GeV, the jet energy obtained at GESR. The three s-states lead to different 
behavior and overall probability. Note that the 1/2 ground state is produced roughly as 
frequently as the 1/2' while the latter is produced about twice as often as the 3/2 state. 
The observed 1/2 ground states are produced from the decays of these vector diquark states 
approximately as often as they are directly produced. 
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In Fig. ^ the corresponding fragmentation functions for charmed states with non-strange 
diquarks are shown. These functions are evolved to 45 GeV, the jet energy attained at LEP. 
Fig. I shows the b-quark states with non-strange diquarks, also evolved to 45 GeV. It is 
clear that these spin-averaged fragmentation functions peak at high z, even after evolution. 
The vector diquark baryon states have noticable secondary peaks at the low scale, which 
get diluted at higher scales.The secondary peaks arise from the polynomials in s in the 
chromoelectric and chromomagnetic contributions and depend on the amount of the latter, 
as hxed by the parameter 1 -|- k. The shapes of the curves are signihcantly different from the 
Peterson function shape and eventually should be distinguishable experimentally. The 
peak position moves towards higher z as the heavy quark mass increases, as expected from 
heavy quark QGD. 


III. SPIN DEPENDENT FRAGMENTATION FUNCTIONS 


In the preceding, the spin orientations of the virtual quark and heavy baryon have been 
summed. In general, however, there is a spin dependence to the fragmentation process. For 
a spin 1/2 baryon there are two fragmentation functions that characterize the leading twist 
spin dependence, gi{z) and hi{z). They correspond, at the parton model level, to the transfer 
from the quark to the baryon of longitudinal polarization (or helicity) and transversity [Q . 
As a consequence of the development leading to Eq. |^, the spin dependent fragmentation 
function gi can be written in the form 

1 lA P_lA P 

gi{z) = DQ^+)^Hi+){z) - Dq(^+)^h{-){z) = Y^^lin^ 

where jA+^P is the probability for the helicity -1-1/2 quark to produce a helicity A baryon, 
with a sum over the unobserved diquark degrees of freedom implied. The Ao+ represents 
the corresponding production of an on-shell helicity -1-1/2 heavy quark, at large go + Qs- 
The reader may ask, how can a heavy quark flip its helicity in hadronization? For the 
scalar diquark combining with the heavy quark, this has to be an effect that would vanish 
in the heavy mass limit, relative to the spin independent fragmentation, i.e. gi{z) would 
coincide with fi{z). For the corresponding vector diquark case this need not be true, since 
the diquark can carry negative helicity leading to the opposite helicity for the baryon. 

The analogous transversity function requires the superposed helicity states (|-1-) ± 

i|-))/V2=|fc). 


hl{z) - DQi^y^)^H(y+)(z) - - 


-- lim 

IGtT"' 90^00 . 


ds 


A. 


y+,y+ 


A. 


y+,y- 


Sth 


A, 


Oy+ 


(27) 


Now in the model we are considering the transversity can not flip, since there is no chirality 
change in the matrix elements. We have the transitions I"*" ^ + O’*" or I"*" -|-1"*", analogous 

to a virtual quark decay into a baryon and a positive parity boson. To get a non-trivial 
result there needs to be an opposite parity bosonic state as well. 

In Figs. ^ and ^ the longitudinal fragmentation function gi{z, Q^) is plotted for the two 
spin 1/2 states with the same physical parameters as in the preceding hgures. Note that 
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gi is very similar in shape to the spin averaged case for the A states, as expected when the 
diquark is a scalar. This shows that the helicity flip contribution is relatively small for the 
values of r = relevent for the bottom and charm baryons. Hence the longitudinal 

polarization of the heavy quark is passed on to the baryon. But we will see that this spin 
preservation is less than 100%. Furthermore, the production of excited baryons will dilute 
that spin preservation. The spin 1/2 S states have very different behavior, showing the 
importance of non-zero helicity for the vector diquarks. 

The spin dependent fragmentation functions for the spin 3/2 baryons are even richer 
in complexity. There are seven such functions at leading twist, many of which will be 
accessible from the decay distributions of these states into the 1/2 state plus a pion. While 
these fragmentation functions have not been classihed in the light-cone expansion formalism, 
it is clear that the number of independent leading twist functions coincides with the number 
of forward amplitudes for parity conserving elastic scattering of spin 1/2 on spin 3/2. In 
general there will be 2(2S-|-1)-1 such leading twist fragmentation functions for a quark to 
fragment into a spin S hadron. In the heavy quark limit the helicity of the quark will be 
preserved in the hadron, so we expect the analog of gi to be near the spin averaged function. 


IV. COMPARISON WITH DATA 

The analysis of charmed baryon production data has been extensive at CESR and, more 
recently, at LEP. Fragmentation data now exist from the CLEO collaboration [Q for some of 
the Sc states. The Ac is studied at both LEP and CESR. Bottom fragmentation into baryons 
is now being studied at LEP. CLEO, in particular, has determined spin independent frag¬ 
mentation functions for the lowest mass spin 1/2 and 3/2 states Sc. Polarization asymmetry 
has been measured for A and Kh at LEP. Production rates for Ac from c-quarks have been 
determined. All of this data provide a testing ground for the model being proposed here. 

The parameters that enter our calculations of fragmentation functions (and their integrals 
over z) are the masses of the constituents, the poles in the chromodynamic form factors, the 
widths of those poles, and the anomalous chromomagnetic moment of the vector diquark. For 
the diquark masses we take m{ud) = 0.6 GeV and m{us) = 0.9 GeV. The baryon masses are 
taken from the data, so the ratios, r = ttid/M oi diquark to baryon masses are determined 
thereby for each baryon. It is assumed that the difference between a baryon mass and the 
constituent heavy quark mass plus the diquark mass is negligible, in order that the Bethe- 
Salpeter wave functions need be evaluated only at the origin. The poles that enter the form 
factors of Eq |T^ for the (ud) diquarks are taken from the electromagnetic form factors 
Ms{ud) = 1.8 GeV and Mv{ud) = 1.2 GeV. The full width at half maximum for the scalar 
form factor is set at 0.88 GeV and, for simplicity, the vector form factor is assumed to have 
the same width. Recall that the width is introduced so that there are no singularities in 
the physical region of s, the virtuality of the fragmenting quark, or, correspondingly, the 
square of the gluon 4-momentum. It will transpire that the fragmentation probabilities will 
depend critically on the pole positions and width, since the integration region is dominated 
by the lowest values of s, where the heavy quark is nearly on mass shell. In that region the 
poles are nearby. The pole and width parameters are expected to be different for the (su) 
diquarks. We take a cue from the electromagnetic form factors of the charged tt’s and iP’s, 
where the charge radii are roughly in a ratio of 1.3:1, corresponding to a pole postion that 
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increases by 1.3 for the strange meson. This is qualitatively understandable by analogy. The 
p vector meson contributes to the tt charge form factor, while the 0 vector meson contributes 
to the kaon electric charge form factor. The latter has a mass 1.3 times that of the p. So, 
to fix the (us) diquark chromodynamic form factors we choose an overall scale factor of 1.4 
for the pole positions and the width, slightly bigger than the meson case. 

The anomalous magnetic parameter hem was determined for the (ud) vector diquarks 
to be a positive number - a result of htting the composite nucleon form factors. For the 
chromomagnetic case, however, a negative value is preferred for k from calculations of the 
mass spectrum of excited baryons Furthermore, using a diquark-quark model of the 
nucleon to calculate the electromagnetic charge radii and polarizabilities preferred a negative 
value for kem- It is unclear what value to take for this parameter, given the divergence of 
different methods. 

We will determine a value for k by optimizing our model predictions compared to data 
for the ratio of production probabilities, R(T,h) = Sb/(Sb -|- S^). In the model, and 

are 6-|-vectord} diquark states of spin 1/2 and 3/2. The difference in production 
probabilities or f\{z) for these two states depends sensitively on n. Using the measured 
value from DELPHI of 0.24 ± 0.12 for the ratio, we choose k = —1.10. This makes the 
overall chromomagnetic coupling small and negative (!-!-«: = —0.10). The reason for this 
small value is that the ratio RiJ^b) would be exactly 1/3 from spin counting if there were no 
chromomagnetic term at all; the 1/3 is compatable with the data. 

The simplest states to study, from our point of view, are the ground states, since they 
involve the scalar diquark. For these states the integral over ^ of fi{z) should correspond 
to the total production probability for producing the state from the corresponding heavy 
quark. However, there are contributions to the same probabilities from the excited states 
that decay into these ground states. Consider the fragmented from a c-quark, for which 
OPAL [^] measures 5.6±2.6% and CLEO p finds 9.5±1.3%. These measurements include 
directly fragmenting Ac’s along with any state that decays into this ground state. The Sc’s 
(both spin 1/2 and 3/2 states) decay strongly into Ac ± tt, so contribute to the rate. With 
the parameters chosen, we hnd 0.5% for the directly fragmented Ac and 3.3% when the Sc’s 
are included. This is consistent with the LEP data. For the analogous b-quark system we 
have hxed the T,b ± rate to be 4.8%, consistent with experiment []^, 4.8 ± 1.6%. The 
total Aft is then 5.8%, comparing nicely with the measurement of 7.6 ±4.2%. These and 
the following results are summarized in the Table below. 

It is signihcant to note that we have obtained these sizeable baryon fragmentation rates, 
in contrast with the similar model of Martynenko and Saleev 0. Saleev obtains only 
0.2% for the Aft. This indicates the importance of our form factors in getting the correct 
normalizations. 

Having conhdence in the overall normalization, we have reason to trust the full frag¬ 
mentation functions. These are shown in Fig. ^ and ^ for the c and b states. The input, 
unevolved “boundary data” show a large peak at high z and a secondary peak at medium z. 
That fairly severe behavior is moderated considerably after evolving to the scale of CESR 
or LEP. But even at the LEP scale, the functions are distinguishable from the Peterson 
function, being peaked at higher and more skewed. As sufficient data is gathered, it will 
be possible to see such a difference. 
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For the singly strange diquark, the lowest charmed baryon 1/2’'' states are the c + [m, s] 
and c+ [d, s] states, and involving the antisymmetric, spin 0 diquarks. These, along 


with the spin 3/2+ states (c + {m, s} and c+ {d, s} baryons), and involving the 
symmetric, spin 1 diquarks, have been seen and measured in sufficient quantities for CLEO 
to sketch their fragmentation functions [^. The spin 1/2+ partners, of 3/2+ states have 
not been seen yet. They are presumed to have a mass below the 5c + tt threshold, so must 
be seen in radiative decay channels. Note that these latter 5^ 1/2+ states have the same 
isospin as the lower lying ground states He 1/2+ and could mix with them, in principle. In 
any case, the measured fragmentation functions provide a crude test of the model. The data 
are £t by the experimenters with a common parameterization of the Peterson function [|^. 
It is easy to see in Fig. that the data fall nicely on fi of our model, evolved to Q = 5.5 
GeV, with the possible exception of the highest z data point. These data are not sufficiently 
accurate to be a crucial test of the model, but do exhibit the trends we expect. Note that 
the experimental variable Xp does not correspond exactly to our z, the light cone variable. 

The ratio of the 3/2 to 1/2 production can be extracted from the data with some uncer- 

:*o ;^+ _|_ 7 ^- jg given 

■ ^ 5° + 7r+ 


tainty [^. The percentage of all 5+ states that arose from decays 
as (27 ± 8)% and the percentage of all 5° states that arose from decays 
is given as (17 ± 6)%. (Note that we have combined the statisitical and systematic errors 
here.) 

The experimenters do not see the channels, H*+ —> 5+ + and 5**^ ^ 5° + 7r°. 
From isospin conservation these channels account for 1/3 of the decays into He + vr, while 
the reported charged vr channels constitute 2/3. Suppose iVH* states of both charges are 
produced. Then 2/3 N will be seen in the charged tt decay mode. The total number of H+’^^’s 
seen will be = §77/(0.27,0.17) (supressing errors until the end). The number of H+’°’s 
not coming from the decays of the 3/2 states will be N^^q — N. Assume that n_|_ o of the H+’°’s 
come from other fragmented states’ decays. Then A^+,o — N — n+o is the number of direct 
fragmentation products of the charmed quark. The ratio i?(+or0) of directly fragmented 
H+’° to S*+’° is given thereby as R{+) = 1.5±0.7—n+/A^: 1 and R{0) = 2.9±lA — no/N : 1. 

The numbers will come from the radiative decays of the heavier 1/2 states, as well 
as higher He states (radial and orbital excitations of the c + (su) and c + (sd) systems). We 
have calculated the fragmentation functions for the spin 1/2 quark-vector-diquark states 
and hence the number of H^ spin 1/2' states vs. H* spin 3/2 states. That is 0.5:1 for the 
parameterization used in Fig. Assuming n+ o is due entirely to these 1/2' states decaying 
100% into the ground state H+’°, we have for the different charge states R{+) = 0.9±0.7 and 
i?(0) = 2.3 ± 1.4, both of which are consistent with the ratio of 1.4:1 predicted by the same 
model calculation. Hence, if the model is taken seriously, and the experimental uncertainties 
are hrm, the data do not require large contributions from fragmentation of the c-quark into 
higher excitations of the He states. 

There are two reasons to be cautious about these experimental numbers, however. First, 
the errors are quite large, leaving considerable variation possible within two standard devi¬ 
ations. Secondly, the CLEO results are obtained at e+ ± e~ energy near 10 GeV. For the 
c-quark jet at roughly 5 GeV, the extraction of asymptotically meaningful fragmentation 
functions is somewhat dubious. 

In a previous version of our model our parameterization gave a much larger vector 
diquark to scalar diquark production probability. With the more reasonable values now 
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adopted these diquark states and the corresponding baryons are produced with roughly the 
same probabilities, as the calculations for heavy-heavy baryons by Martynenko and Saleev 
favored. 

Finally we consider the spin dependent fragmentation. At this time there is not enough 
data to determine 2 ; dependence for polarization in heavy quark fragmentation. However, 
there is a determination of the net longitudinal polarization of A;, produced at LEP. That 
number is —0.23 ± 0.25 as determined by ALEPH |^, using a technique suggested by 
Bonvicini and Randall Given that the b-quark produced at the pole is expected 


to have longitudinal polarization of —0.94, this measurement gives 0.24 ± 0.27 for the net 
transfer of helicity from the b-quark to the A;,. This is rather low if one anticipates that 
the heavy baryon carries most of the helicity of the heavy quark - the expectation of heavy 
quark held theory . 

Now the longitudinal polarization of the A^ fragmenting from a positive helicity b-quark 
is a function oi z - the ratio gi{z)/fi{z) for the directly produced A;, (calculated to leading 
twist). The integral over z of that ratio would give a net polarization. Experimentally, 
though, the net polarization is obtained by taking each event, regardless of its 2 ; (and pr), 
and calculating a quantity related to its polarization. The result of this process is to give 
a net polarization that will be the integral of gi{z) over divided by the corresponding 
integral of fi{z). Both of these integrals are scale independent. Using our fragmentation 
functions we obtained 0.90, only marginally lower than the heavy quark limiting value, but 
still not the small result extracted from the data. However, from our spin independent 
calculation above and the LEP data, we know that the S;, and are produced in relative 
abundance, and will decay into A;, -|- vr, so that the polarization will be diluted by these 
other channels. A heavy quark limit calculation by Falk and Peskin ^2\ anticipated this 
circumstance. They summed the contributions of all of these states in determining the net 
Afo polarization. Using simplifying assumptions, they obtained about 0.72 for the fractional 
helicity transfer from the b-quark to the hadron. That estimate assumed the ratio (called 
A in their paper) of Sb and (vector diquark states) to A;, (scalar diquark states) of 0.45. 
Taking that ratio to be 4.8 instead, which is our result and consistent with experiment, the 
resulting fractional helicity transfer becomes 0.26. Falk and Peskin also dehne a parameter 
Wi which measures the amount of helicity ±1 diquark that combines with the heavy quark. 
They take that parameter to be zero, whereas we can calculate wi, which averages over all 
.2 to be ~ 0.6. With our values for both A and wi we obtain the fractional helicity transfer 
of 0.46. Both of these results are near the central value of the measurement. In obtaining 
these results we have assumed the heavy quark limiting values of ±1/3 for the polarization 
of the secondary A^’s resulting from the Sb and decays. These results will obtain when 
the chromomagnetic contribution is small, as it is for the k we have chosen. 

The application of a similar model by Saleev |l^ to Ab production and polarization 
yields very different results. As we noted above, his production probability is too small. 
The polarization for direct Ab is less than ours (0.6 to 0.7). However, it seems that it is the 
chirality asymmetry (Left - Right) that Saleev has calculated, rather than helicity. It is the 
latter that is measured experimentally. We have a much weaker dependence on the diquark 
mass in our calculation as a result. 
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V. SUMMARY 


A model for fragmentation into heavy flavored baryons has been developed using pertur¬ 
bative QCD and a Bethe-Salpeter wave function for a quark-diquark system. This provides 
the starting point for QCD evolved fragmentation functions. The parameterization of the 
diquark structure through the chromodynamic form factors for scalar and vector diquarks is 
accomplished by using the pole form applied to the electrodynamic form factors. The kine¬ 
matic region near these poles is quite important because the integral over the fragmenting 
quark’s virtuality, s (or indirectly, the baryon’s transverse momentum at a hxed z), empha¬ 
sizes the nearly on-shell region where the corresponding gluon is near the poles. It is this 
pole parameterization that is crucial for determining the overall magnitudes of the various 
production probabilities. 

The production probabilities were all close to experimental values, which supports our 
reasoning about the form factors. The and dependences predicted have the common 
feature of being very sharply peaked at high 2 : for the input scale juo, and more broadly 
peaked at high Q^. For the 1/2+ ground states, the dependence of the spin dependent gi 
is close to the form for /i, but their ratio (which will determine the baryon polarization as a 
function of z) is striking. For the higher mass 1/2''' state containing the vector diquark the 
gi has a very different 2 : dependence. It will be particularly interesting to see if the peak 
and dip structures for both spin 1/2 states are reproduced by the data. 
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Particle 


P{c 

Ac) (including decays of Sc and S*) 

P{b 

^ Afe) (including decays of S^ and S^) 

P{b 

-s, + s,*) 


P{b^^b) 

TJJP 

>S5)+P(6^S:) 

P(c- 




P{c 

^ He) (including decays of H^, and H*) 

P{c 

^ He) (direct production) 

P{c 

^ H*^ or H**^) (direct production) 

P{c 

^ Ac) (direct production) 

P{c 

^ Sc++ , S*+ or S*°) (direct production) 

P{b 

^ Afe) (direct production) 

P{b 

^ S^^ , or S^“) (direct production) 


Prediction 

3.26% 

5.8% 

4.8% (fixed) 

0.33 

1.4 

0.53% 

0.17% 

0 . 12 % 

0.52% 

0.62% 

1 % 

1.06% 



































